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1 Introduction 



With the new colhder facihties, upcoming experiments will bring higher statistics 
data samples into the low energy regime. Due to their accuracy higher order 
calculations are needed to update the theoretical prediction for the measurements. 
In this frame the calculation of the two-point functions at next-to-next-to-leading 
order (NNLO) at low energies has become necessary. These provide us with the 
pseudoscalar masses and the decay constants, which are needed input for most 
other quantities. The two-point Green functions are also basic tools in the study 
of the strong interaction. They form the basis for a series of very useful sum 
rules starting with the Weinberg and DM0 sum rules (we refer to [Q and 
references therein for a more complete discussion about sum rules). 

In this work we are concerned with the low energy regime of QCD. We will 
study the two-point functions with Chiral Perturbation Theory (CHPT), valid for 
energies below the first resonance {E ~ mp) and describing the strong interactions 
using the pseudoscalar octet as the basic fields. This is by now a fairly developed 
field. We refer to for reviews and various abstracts on recent works. 

For a future study of various sum rules, we present the vector and axial- vector 
two-point functions at NNLO in three flavour CHPT, in the limit of unbroken 
isospin. Four of the six basic correlators have been calculated earlier |^, |^ and 
we have fully conflrmed their results in the vector-vector two point function and 
partially[| in the axial-vector case. The other two are new and complete the three 
flavour basis. As a byproduct we also give the masses and decay constants to 
NNLO. 

The interest in the NNLO calculation is beyond the precise measurement of 
the couplings and masses of the effective theory, it allows to test the convergence 
of the theory and provides a more stringent check on the principle of resonance 
saturation of the constants in the low energy chiral Lagrangian. While this prin- 
ciple worked well at 0{p'^) ^ only a few tests at 0{p^) have been done. In 
this paper we estimate some of the constants appearing to two loops and check 
their effect on the full 0{p^) results including the loop contributions. We use the 
MS-subtraction scheme and the recent classiflcation for the 0{p^) Lagrangian ||^. 

Some applications to chiral sum rules for the isospin and hypercharge cases 
can be found in 0, [l^]- We intend to return to that subject in a future pubhcation. 

The paper is organized as follows. In Sect. || the two-point Green functions are 
deflned, followed by a short overview of CHPT in Sect. |^. The vector two-point 
function is relatively easy since it only involves products of one-loop integrals. 
Its calculation is described and results given in Sect. ^ The masses and decay 
constants can be calculated in two ways. The masses can be obtained from 
the zero of the inverse pseudoscalar propagator or from the pole in the axial- 

^ The authors of Refs. ^ use a different method to perform the sunset integrals making 
comparison of those parts difficult. We agree on all the parts we could check without converting 
their sunset functions to ours. 
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vector two-point function. The decay constants can be directly determined from 
the residue of the poles in the axial-vector two-point function or through the 
definition with the axial-vector-pseudoscalar function. We have checked that 
both methods agree and the first is described in Sect. |. Finally the axial- vector 
two-point functions are presented in Sect. |. The new OijP) constants appearing 
are estimated in Sect. ^ on the basis of Resonance Dominance. In Sect. § some 
results are presented, postponing a more detailed analysis [|ll|]. And finally in 
Sect. 1^ we discuss our main results. 

We refer the lengthiest expressions and the more technical discussion of the 
loop integrals and renormalization to appendices. In App. |A.1| and App. [A. 2 



we 



give the full expressions for the masses and the decay constants. In App. P we 
display the the axial-vector two point function components. The loop integrals 
are collected in App. 0. 



2 Definition of Two-point Functions 

We calculate in CHPT the two-point functions of vector and axial- vector currents. 
The quark currents are defined by 

V;Kx)=il^.(l' and 4^(a;) = r7^75g^ (1) 

where the indices i and j run over the three light quark flavours, ti, d and s. 
Working in the isospin limit all SU(3) currents can be constructed using isospin 
relations from 

V;{x) = i=(\/"(a:)-Ff(x)) , 



y,{^) = V-^\x). (2) 

We refer to these as isospin, hypercharge and kaon respectively. E.g. the electro- 
magnetic current corresponds to 

We calculate the two-point functions defined as 

^^Jd'x e^«- (0|T(V;(x)Kr(0))t|0) , (4) 

for a = 71,7], K. All other vector two-point functions can be constructed from 
these using isospin relations. Lorentz-invariance allows to express them in a 
transverse, and a longitudinal, II'^''^ part 

= - 9'^7..)nK(g^) + g,g.nS.°l(g2) . (5) 
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Similar definitions and comments apply for the axial- vector currents. 
The currents obey Ward identities and other symmetry relations; for = 
the vector currents are conserved 

T^P.iQ') = = . (6) 

In the SU{3)v limit, the three vector two-point functions reduce to the same 
expression. The same holds for the axial-vector currents. In the addition the last 
ones are conserved only when the relevant quark masses vanish. 



3 Chiral Perturbation Theory 

Effective theories are a general tool in understanding a wide range of physical 
processes, from high energy physics to superconductivity. In that frame Chiral 
Perturbation Theory is a successful theory describing the strong interaction at 
low energy. It is based on the existence of a mass-gap in the hadronic spectrum, 
at low energies only the low mass states can be excited. Those are the Goldstone 
boson particles and are the only states that are actually predicted from first 
principles in QCD [|l^]. For the present status of the field we refer to the listing 
of review articles and the abstracts in and to some recent lectures . 



For constructing the effective action, the high energy states of the theory 
should be integrated out, thus the Lagrangian describing processes at low en- 
ergy consists of a series of operators involving only Goldstone boson particles. 
These general operators should share the same symmetries as the basic underly- 
ing theory [|l^, |15|, in particular Lorentz invariance, local chiral symmetry, parity 



and charge conjugation. Thereby in the following we use that the generating 



functional of both theories, QCD and CHPT, are the same at low energies |15 



exp{iZ} = j VqVqVG^ exp |z j d^xCqcD^ = J exp |i j d'^xCefj^ ■ (7) 

Following this philosophy, the QCD effective Lagrangian is given by a non-linear 
realization of chiral symmetry (see [|l^ and references therein). The lowest order 
in an expansion by quark masses and external momenta is 

A = ^{D^UD^U^ + U\ + x^f/), (8) 

where 

f/(0) = = exp(2y2$/Fo) , (9) 

parametrizes the pseudo- Goldstone bosons and 

X = 2Bo{s + tp) (10) 
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is given in terms of the scalar and pseudoscalar sources s and p. {X) = tr flavour {X) 
and U,s,p are matrices in flavour space. Both Bq and Fq are constants not re- 
stricted by symmetry. They are related with the quark condensate and the meson 
decay constant respectively. To respect local invariance the external sources are 
incorporated through the covariant derivatives 

D^U = d^U -ir^U + iUl^, (11) 

and the field strength tensors 

Fl" = d^F - d'^F - 1% F^^ = d^r'' - d-^r" - z[r^, r'^] . (12) 

As mentioned above the purpose of this paper is to compute two-point Green 
functions, the most straightforward way is to incorporate classical sources in the 
effective action. This reduces the calculation of any n-point Green function to 
the evaluation of functional derivatives acting on the generating functional 



G^-\x^, ...,x^)= Z[J] 

6j{Xi) . ..6j{Xn) 



(13) 

J=0 



This allows for instance, to relate the chiral condensate with the constant Bq by 
taking the derivative respect to the scalar sources. This formalism is not only 
suitable for an easy calculation but also allows to incorporate the electromagnetic, 
weak interactions and the symmetry breaking through the quark masses via the 
following identifications 

r^^ = Vf, + af„ l^ = v^- a^, s = M + . . . , (14) 

where M. stands for the diagonal quark mass matrix, M. = diag {rriu, rrid, rris) . 

To get the desired chiral order — 0{p^) — in our calculation we will deal with 
three kind of contributions: i) tree, one and two-loop diagrams involving vertices 
of C2, a) tree and 1-loop graphs with vertices from C2 and from £4 given by 

^i(V^f/^V^f/)2 + L2(V^f/1-V.f/)(V^?7tV"f/) 

+L,{V,UW^U{x^U + xU^)) + L,{x^U + xU^? 
+Lr{x^U - xU^? + L,{x^Ux^U + xU'xU^) 
-iL,{F^^V^UVM^ + FrV^U^VM) + L.oiU^F^'UFL,,) 
+Hr{F^'FR,, + FI'Fl,.) + //^(x^x) (15) 



and finally Hi) the tree graphs of Cq. The latter was first classified in [17] and 
recently in |]^ a more restrictive general set was found. We borrow in the following 
our notation from the last reference. Notice that the terms Hi and H2 in Eq. 
(|T3p have no direct physical meaning, their value depends on the precise way in 
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which the currents are defined in QCD. But once a consistent definition of a QCD 
current has been given they are defined unambiguously. 
For later use we define the following quantities 

= (Fr + Fr)/2, 

= ^{u\d^-irf,)u + u{df,-ilf,)u^} , 

X+ = u^XU^ + ux^u , 
V^X = d^X + T^X-XT^. (16) 



4 The Vector Two-Point Functions 

Within the framework of previous sections we can start to calculate the vector- 
vector two-point functions. The first contribution appears at 0{p^). It has been 



calculated for the isospin case in [T^ in SU(2) CHPT. The extension to 0{p ) 
in SU(3) CHPT has been done in for the isospin and the hypercharge case. 
We have reproduced their results and in addition we present the kaon vector 
two-point function as well here. 

The calculation to 0{p^) presents no new difficulties besides being rather 
tedious since only products of one-loop integrals appear. 

The O(p^) contributions come from diagrams (a-c) in Fig. |l|. Diagrams 
(d,e,m) and (n) can be calculated directly or using wave-function renormalization 
and mass corrections. We have checked that both approaches give the same result. 
As a consequence the final result contains no three-propagator integrals. They 
cancel when the O(p^) result is expressed in terms of the physical masses using 

/K; g^) = fiml, q^) + - ^?)ct/K; ^')L. (17) 

where the are the bare masses and the mf the next-to-leading order masses. 
In addition we replace Fq by and all masses by their physical ones in the 0{p^) 
expression. 

There are no one-particle reducible contributions to the vector two-point func- 
tions. 

We have performed the following checks 

1. In the isospin and hypercharge case the longitudinal part vanishes. 

2. In the SU(3) limit, i.e. = = ^Ihirah two-point functions 
are equal. 
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(d) (e) (f) 




(g) (h) (i) 




0) (k) (1) 




(m) (n) (o) 

Figure 1: Diagrams for the vector two-point function. The crossed circle stands 
for the 0{p'^) vertex insertion. Wiggly hnes are the external vector currents. 
Dots are C(p^) vertices and a square is a 0{p^) vertex. The solid hnes are meson 
propagators. 



6 



3. The SU(3) breaking effect in the form-factors appears only in second order 
in the quark masses, i.e. order (m^ — m)^, as required by the Ademollo- 



Gatto theorem 19 



4. All divergences with a non-analytical dependence on masses or cancel 
and the ln(47r) and 7 terms can be absorbed in the counter-terms as well. 
Both of these follow from general renormalizability theorems. 

5. The remaining divergences are in agreement with those from the general 
calculation using heat-kernel methods and with the double logarithms 
of ill. 

The result can be expressed simply in terms of the finite functions 



for 2,j = 71, K,!]. The definitions of B22{Tnf,m'j,q'^) and Aq can be found in 
App. |C1| . Notice that By{q^) is regular at = for i = j. 
For the isospin transverse part we find 

n« = -4(2/7[ + L[,) - 8i?r(g') - 45^ ^(g^) 

(19) 



The hypercharge transverse part is given by 



n« (g^) = -4(2iJ[ + L\,) - 125^^(g2) 



- ASq^B^^{q')Ll + 48/ix(i^S + ^lo) + ^Qq\B^^ {q')f^ 
-^2mlCl, - 32(m2 + 2ml)Cl, - Sq^C;, . (20) 

The longitudinal part vanishes for the above two. These results agree with those 
obtained in |^ when the differences in subtraction schemes are taken into account. 

The expressions for the kaon two-point functions are new and are somewhat 
longer. The transverse part is given by 

4k?') = -mHl + Ll,) - 6{B-y^{q') + Bf{q')) 
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1 r 8 



4g2 



32 



-32m|Ce^, - 32(m2 + 2ml)Cl, - ^(m^ - mlfC;, - 8q'C;, . (21) 



This two-point function has also a longitudinal part 



nVl(g^) = ^((m^ - mirr^iq') + {ml - m\fj-^{q^ 

+ ^{ {^{rn\ - ml){3fi^ - 2fiK - f^,) + \2{m\ - mlfr^{q^) 
_ \ 3 

3 

+l2nKHn - ^fJ'rifJ'ri) - ^(-5/i,r + '^fJ'K + 3/i^)(m|- - ml)r^{q^) 
9 / ^2 



16 

+ 



« - miyr^iq^) + {ml - m^fj^^iq^) 



16g 



+-[mi-ml)\ACl, + Cl,), (22) 

Notice that the Ademollo-Gatto theorem is explicitly satisfied. 

All divergences have been absorbed in the coefficients of the (9(p^) Lagrangian 
by setting 



4 3 

/X ^Cgi = - ^^Ai(L9 + L 
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(a) (b) (c) 



(d) 



Figure 2: The diagrams contributing to the axial-vector two-point function. 
The filled circles indicate the one-particle-irreducible diagrams. Solid lines are 
pseudoscalar meson propagators and the wiggly lines indicate insertions of an 
axial- vector current. For the inverse propagator the wiggly lines are meson legs 
and for the decay constant the right wiggly line is a meson leg while the left 
remains an axial current. 



1 , 1 



(1 - C93 - ^^^-^2 + ^^-^1-^9 ' (23) 



which agrees with the calculation of . 



5 Masses and Decay Constants 
5.1 Masses 

The definition of mass is the position of the pole in a two-point Green function 
G(p^,m^) that contains the relevant particle as a possible intermediate state. 
The axial-vector two-point function is a suitable candidate to obtain the masses 
for the pseudoscalar mesons. The general structure of this two-point function 
is shown in terms of one-particle-irreducible (IPI) diagrams in Fig. |^. However 
the meson-propagator or the two-point function of meson fields itself is a simpler 
set with the same pole. We denote the sum of IPI graphs by in(p^,m^Q,F( 



0) 



The set of diagrams contributing to n(p^,m^Q,Fo) is depicted in Fig. ^(c). The 
propagator is 



\ n 



-ml + Il{p^,m%,Fo) 



(24) 
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(a) (b) (c) 



.0. 



(d) (e) (f) 



.0- 



(g) (h) (i) 



.0. 




(j) (k) (1) 

Figure 3: The set of diagrams contributing to the IPI quantities, (a) axial- 
vector-axial- vector (b) axial-vector-pseudoscalar (c) pseudoscalar-pseudoscalar. 
(d)-(l) the respective diagrams when the dashed lines are replaced with the exter- 
nal legs of (a), (b) or (c). A line is a meson propagator, a wiggly line an external 
source, a dot a vertex of C(p^), a crossed circle a vertex of 0{p'^) and a crossed 
box a vertex of 0{p^). 
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where ml stands for the lowest order mass and m^Q collectively denotes the various 
lowest order masses. The physical mass is given by the zero of the denominator 
once the external legs are on mass-shell 

^Ihys -^1 + ^{^Ihys, "^L ^O) = . (25) 

We replace the masses rrijo by their physical masses and Fq by F.^. It is sufficient 
to use the NLO formulae for these in Il{p'^,m'^Q, Fq) since it is already of O(p^). 
This leads to 

m%y, = ml + (mT^ + {m'Y'\ (26) 

where the bare masses appear only at the leading order and superscripts refer to 
the chiral order. 

The resulting formulae for the pion, kaon and eta masses are gathered in App. 



A.l| . The formulae for the pion and eta mass agree with those of (the explicit 



formulae only appear in the preprint version) in the way described earlier, see 
footnote 0, while the kaon result is new. 

Notice that the precise expression for the 0{p^) is dependent on the choice of 
the 0{p'^) expression]^, using the Gell-Mann-Okubo relation at C(p^), produces 
differences at 0{p^). 

The masses depend on seven combinations of the 0{p^) constants. All the 
relevant checks described in Sect. ^ were done. 



5.2 Decay Constants 

The pseudoscalar decay constants. Fa, are defined by 

{0\A^{0),\r{p))=^V2p,Fa. (27) 

They can be obtained directly from the definition or by the residue of the pole 
in the axial-vector two-point function. We have calculated them first using their 
definition and verified that the calculation via the two-point function yields the 
same results. 

This calculation involves the use of the expression for n(p^, mf^, Fq) obtained 
earlier to get the wavefunction renormalization in addition to those diagrams of 
Fig. ^(b) for the matrix element itself. 

We can then write the results in the form 

Fa = Fq{i + Fi'^ + )) . (28) 

Similarly to the masses, the precise form of Fj^^ depends on the choice of F^^\ 
For and we have checked the double logarithms with those presented in 
and the result for Fx is new. 

The explicit formulae are rather long and can be found in App. |A.2| . The 
relevant checks described in Sect. ^ were performed. 

^ This is why our expressions have some differences with those of |^ even after correcting 
for the renormaUzation scheme. 
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6 The Axial- Vector Two-Point Functions 



The axial- vector two-point functions to lowest order are quite simple and all three 
reduce to 

2_p2 

n^."(g) = 2f2^,. - -r^w^ ■ (29) 

The NLO corrections only introduce minor changes. The decay constants change 
to the masses to the physical ones and there is an additional contribution 
from the 0{p'^) constants to n^^(g^) = 4L^q — 8if[. In fact a very large part of 
the 0{jP) corrections is of a similar nature. We thus define 

n;^,"(g) = 2Flg,, - -^^q^q, + {q^q, - g^,q^)Il^}l{q^) + q^q^IlfM'') ■ (30) 
q 

The function n^^(g^) can be fully calculated from diagram (a) of Fig. ^. These 
are depicted in more detail in Fig. ^(d-l) and discussed in App. [B.l , 

All the diagrams in Fig. |^ contribute to n^^(g^) even though most of their 
contents actually go into the redefinitions of the respective decay constants and 
masses. The full result is given in App. p3.2| . The results fulfill the same checks 
as in Sect. ^ We call n^](g^) and n^^(g^) the longitudinal and transverse 
remainder respectively. 



7 Estimates of some 0{p^) constants 

In this section we estimate some of the 0{p^) constants that appear in the results. 
We assume saturation by the lightest vector, axial-vector and scalar mesons, 
extending the formalism used in to the present case. 

For the spin-1 mesons we use the realization where the vector contribution to 
the chiral Lagrangian starts at 0{p^). Keeping only the relevant terms for our 
calculation we have 

>Cv = -\{V,,V^'' - 2M'yV,Vn - ^( V/D + • • • , (31) 

where 

and the same holds for the axial- vector with the label change V ^ A and /+ ^ 
/_. and are three- by-three matrices in flavour space and describe the full 
vector and axial-vector nonets, thus we assume a nonet symmetry throughout 
the rest of this section. The rest of the notations was already presented in Sect. 

I- 
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After integrating out the vectors the terms contributing to the two-point 
Green functions at 0{jP) are 



^^VA = 4^(/rV,V./f ) + ^^{rV.VJ^^) . (32) 



In the scalar case, the Lagrangian reads 

Cs = ^{WSV,S-MlS') + c,{Su^u,) + c^{Sx+) + ^{S\+) 

+c,{su,j'f) + c:,{sf.,jn ■ (33) 

After integrating out the scalars, the 0{p^) contribution we are interested in 
comes from the terms 

^5 = ^(X+V + + + (34) 

obtained after the shift of the vacuum expectation value and using the equation 
of motion for the scalars. Note that only the relevant terms are written and, as 
in the vector and axial full nonet of scalars is assumed in S. 

As input parameters we use 

Mv = = 0.77 GeV, /y = 0.20, = = 1.23 GeV, = 0.10, 

Ms = 0.93 GeV, = 0.042 GeV, q = 0.032 GeV, 

= 19 ■ 10"^ GeV-\ c'^ ~ and = -2.4 . (35) 

Ms and dm are obtained from the masses of the scalars i^'Q(1430) and ao(980). 
The value is obtained from r(/o — >■ 77) = 0.56 ± 0.11 keV and is compatible 
also with r(ao — >■ 77) = 0.30 ± 0.10 keV. The values of Cm and Cd are obtained 
forcing the saturation of some of the O(p^) constants by the scalars [0 and are 
compatible with those obtained in T{aQ — > tti]). c'^ value cannot be determined 
from data at present, we assume a value similar to c^. 

Using the notation of for the 0{p^) terms, the spin-1 Lagrangian yields 



C's.~l(^-^J-7-e-lO-'GeV-^, (36) 



f2 

Co. ~ -17 ■ 10-^ GeV"' , (37) 

93 ' V ; 



and the scalar Lagrangian estimates 



^ ^ _4.3 ■ 10-3 GeV-2 , (3^ 

Mg 
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^38~TT774 ~l-2-10"' GeV-^ (40) 



2M^ 

Ck-'-^-l.9-l0-' GeV-^ (41) 



Qo~^~l-9-10-^ GeV-^ (42) 

We stress that the goal of this section is to roughly estimate the values of 
the 0{p^) constants. A real determination would imply the use of chiral sum 
rules or other processes to fix them. It is worth mentioning that the result in Eq. 
(|37D is the same if we use an antisymmetric formalism for the vector Lagrangian, 
and is in agreement with the result extracted from the experimental data using 
sum rules for the vector- vector two-point functions ^ . The generalization to the 
three flavours introduces a new relation of 0{p^) constants, 4:C^^ + Cg^, due to 
the explicit chiral symmetry breaking in the kaon Green function. 

Finally, we remark that the precise value of the constants can have an impor- 
tant variation depending of the input values in Eq. (pS]). Consequently, although 
the values cited in this section are used for the numerical results, with the under- 
standing that the other counter-terms are set to zero, we have to keep in mind 
that these values could overestimate the physical ones. The latter is especially 
true for dm since the Ko and mass difference appears unnaturally large. 



8 Some numerical results 



We defer a more accurate comparison with experimental data to [jrT|, but we 
would like to present some results using our explicit expressions. We use the 
values for the Cf obtained in the previous section and two sets of the C(p^) 
constants. They only differ in the values of L\, Lg. Set A is obtained 

from the fit of the unitarized calculation while set B refers to and vrvr 



data at one loop accuracy ||2^. We give both sets to show an example of the 
variation with the 0{p^) constants. We do not show results for varying the other 
LI but this results in a similar variation in size of the 0{p^) results. The explicit 
values we use, at /i = 0.77 GeV , are 



LI 


= 0.37•10~^ 


Ll = 1.35- 10"^ 


J-r _ 


-3.5 • 10~3 


set A 


L\ 


= 0.60 ■ lo-^ 


= 1.50- 10-^ 


Ll = 


-3.3 • 10-^ 


set B 


LI 


= -o.3■lo-^ 


LI = 1.4- 10-^ 


l; = 


-0.2- 10"^ 




jr 


= -o.4■lo-^ 


LI = 0.9- 10-^ 




6.9 • 10"^ 




-^10 


= mfv-fi) 


= -7.5- 10"^ 








H[ 


= mfv + fi) 


= 6.25- 10"l 









(43) 
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Figure 4: Different contributions to the real part of the vector- vector two-point 
function in the isospin and the hypercharge case. The label CT indicates the 
contribution from only the counter-terms L[ at 0{p^) and C[ at 0{p^). jy" gives 
the C(p^) result. VMD indicates the vector model of Eq. (^). 

For which can not be obtained experimentally, we take the value from the 
Meson Saturation Model. Because the vector contribution should cancel for the 
axial- vector two-point function, we use the same model value for LJq. 
The rest of the quantities we use are 

= 0.135 GeV, tuk = 0.495 GeV, m^^ = 0.548 GeV and = 0.0924 GeV. 

(44) 

They seem reasonable averages of the various isospin related ones. 
8.1 The Vector Two-Point Functions 

In Fig. ^ and Fig. ^ we plot the real part of the three vector-vector two-point 
functions choosing set A inputs. 

In all the three cases the slopes are given mainly by the 0{p^) constants 
as estimated above. We have shown these contributions in the curves labeled 
(CT). Essentially — and with exception of C93 — the main effect of varying the 
input parameters is to shift the plots vertically. We see that the loop effects are 
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Figure 5: Different contributions to the real part of the kaon vector-vector two- 
point function. The superscripts L (T) refers to the longitudinal (transverse) 
part. For comparison we also plot the isospin case. (CT) as in Fig. |[ 

larger in the isospin case and smaller for both, the hypercharge and kaon. In 
the chiral limit all three cases reduce to the same, and differences are related 
to the breaking of the SU{3) symmetry. For the isospin, the two-pion channel 
produces the notable difference with the counter-term contribution, while for the 
hypercharge and kaon the smaller difference is explained by the larger masses in 
the loops and some explicit breaking of the symmetry through the quark masses 
in the counter-term contributions. 

In Fig. ^ we also plotted the case with a complete saturation by the vector 
meson (VMD) 

M^v)--—^r-^, (45) 

with r = 0.150 GeV. The conclusion is that models with only vectors explain the 
main part of the two-point function, however an important contribution coming 
from the two-pion intermediate states is present. The curve including only the 
counter terms — isospin (CT) — can also be obtained with the first two terms 
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set A set A set A 



set B 



H (GeV) 

Fk/F^ 
(Fk/F^) 
Ffj/ F.j^ 

phy 



ml /ml 



m 



K 



mi/m^ 



rjphy 



0.77 
0.068 

0.216 

0.312 
-0.039 
-0.003 
-0.045 



0.77 
-0.101 
(0.013) 
0.055 
(0.08) 
0.092 
0.214 
0.241 
0.312 



0.5 1.0 
-0.066 -0.100 

-0.023 0.100 



-0.011 
0.132 
0.246 
0.234 



0.150 
0.238 
0.194 
0.273 



0.77 
-0.172 
-0.050) 
0.035 
(0.06) 
0.065 
0.355 
0.423 
0.521 



-0.001 

-0.450 

-0.600 
-0.003 
-0.873 
-2.428 



Table 1: 0{p*) and 0{p^) contributions to the masses and decay constants. 



The number in brackets are the extended double log approximation of Ref. [pi 
Columns three to six are the 0{p^) loop contributions only. The contribution 
from the C[ is listed separately in the last column. 



of the expansion in the previous formula with F = and considering the tiny 
modification due to the scalars. 



8.2 Masses and Decay Constants 

We continue our discussion with the masses and decay constants. We have sum- 
marized our numerical results in Table |I| using the values for the O(p^) constants 
quoted above. As one sees in columns three to six, both masses and decay con- 
stants have substantial 0{p^) loop contributions. In addition the pure polynomial 
piece at 0{p^) tends to have the opposite sign and is very large using our model 
dependent estimates. This only reinforces the statement in Sect, ^of the lack of 
knowledge in the scalar sector. The terms containing are the only ones con- 
tributing in this subsection and seem severely overestimated even though they 
are of a size expected by naive dimensional analysis. 

In order to have a full presentation of the 0{p^) contributions a refit of all 
0{p'^) coefficients using the full 0{p^) expressions would be needed. We postpone 
this till after the main other processes are also calculated to this order given the 
dependence of the 0{p^) contributions on L\-U^. As an example using set A at 
/i = 0.77 GeV otherwise but shifting to 1.0 ■ 10"'^ reproduces the experimental 
value of Fx I Ft, when setting (i^ = 0. 

For the decay constants the (9(p^) contributions to the ratios are smaller 
than the 0{jy'\ not including the estimates from scalar exchange to the (9(j9^) 
constants. 

To judge the effect of the (9(p^) contributions on determining the quark mass 
ratios we use the lowest order, O(p^) and (9(p^) formulae in terms of physical 
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quantities to obtain the lowest order masses using Eq. p6|). This leads to 



2mo 



rriu + TUd 

2ms 
rriu + rud 



Sml^ — m, 



2 

Ott 



25.9(p^); 24.9(/); 24.1(/) and 



24.2(/); 24.4(/); 23.3(/ 



(46) 



using the results from set A at /i = 0.77 GeV and C, 
compared with 2ms/(m„ + rrid) = 25.5 EH]; 22.8 [El 



= 0. These ratios can be 
extracted from QCD sum 
rules and lattice calculations respectively. 

The emerging conclusions about the convergence of the chiral series should 
be very cautious since a full study includes also the effect of the 0{p^) constants. 
However, while the corrections calculated are significant they do not show evi- 
dence of a breakdown of the chiral expansion for the quantities presented here. 



8.3 The Axial- Vector Two-Point Functions 

In Fig. P we plotted the dependence on momenta of the real part for the re- 
mainders of the axial-vector two-point functions for the three cases under study. 
Because both the longitudinal and transverse remainders have poles at = 0, we 
show the combination fl^'* + fl^^ . A priori we would expect a different behaviour 
for the isospin case due to the three pion channel. However there is virtually 
no effect because the imaginary part is very small in the energy region we are 
considering in agreement with the dominance of the ai axial meson. For the 
other two cases even the three pseudoscalar channel is far. The curves are thus 
very linear. The vertical shifts are due to the explicit breaking from the quark 
masses. The 0{p^) contributions are rather small, the scale in the plot should be 
compared with |2F^/g^| which is larger than 0.07 for the entire region plotted. 



9 Summary and Conclusions 

In this paper we have calculated to NNLO in CHPT the vector and axial-vector 
two-point functions in the isospin limit and in the complete three flavour basis. 

In the vector-vector case, we confirm previous results for the isospin and 
hypercharge 0. 

For the axial-vector case, besides the cancellation of the non-analytic poles, 
we obtain the same double and simple poles that appear with the use of the heat 



kernel expansion pO[. We also agree with the double logarithms, appearing in 
previous work , for the isospin and hypercharge cases. All these checks give us 
confidence about our result. 

We have also given expressions to NNLO for the masses and decay constants. 
The Lagrangian at 0{p^) contains a rather large number of free constants. We 
have estimated some of them using a simple resonance estimate and used this 
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< 



0.025 
0.02 

0.015 
0.01 

0.005 


-0.005 
-0.01 

-0.015 
-0.02 

-0.025 



-0.03 



T 1 1 r 



lsospin<'^-''-'(CT) 



Kaon<"^*'-*(CT) 



Hypercharge<'^'''-'(CT) 
lsospin<'^+'-' 



Kaon 



(T+L) 



Hypercharge*^"^'"' 



-0.25 -0.2 -0.15 -0.1 -0.05 0.05 0.1 0.15 0.2 0.25 

(GeV^) 



Figure 6: Different contributions to the real part of the remainder axial-axial 
two-point function. The superscripts L (T) refers to the remainder longitudinal 
(transverse) part. (CT) as in Fig|^. 

to present some first numerical results for the two-point functions, masses and 
decay constants. We also studied somewhat the /i-dependence of the final result. 

Although the corrections are significant they do not show evidence of a break- 
down of the chiral expansion. For instance, our estimates of the quark mass ratios 
are in agreement with previous determinations. However, the sensitivity to the 
input values, indicate that the C(p^) constants need to be refitted using the full 
0{jP) expressions and that better estimates of the 0{jP) constants are necessary. 
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A Explicit results for the masses and decay con- 
stants 



A.l Masses 



The masses are split as follows 



loops ■ 



(47) 



with moa the contribution from the bare masses. In the 0{p^) we have explicitly 
separated the chiral loop contribution from the model dependent counter-terms. 

In a previous step the masses are obtained in terms of only the bare masses 
(quark masses), we rewrite the 0{p^) contribution with the physical masses im- 
plying a modification of the 0{p^) terms. In the 0{jP) we can safely replace bare 
masses with physical masses. 

For the pion we obtain 



+Umlm\Cl^ - Umlm],C[, + Umlm],Cl^ + 192m>^C2^, , (49) 



+ (56/27m^ — 16/27m^ + 86/27m^m^ ) LU - (589/36 ml + 2 m%) fil 



-Amlfi^fiK + 5/3 m^/i^/i^ - 112 m^/i^L^ - 64 m^/i^L^ 

-56 ml /i^ LI + (144 m^ + 80 ml m\) /x, + 96 ml /i^ LI 

+ (-256 ml - 160 ml ml)fi^ L^ - 176 m^/i^ LI + (-2/3 ml m^^ _ 7^2 ^ ^2^ 

-A/3mlnK /ir, - 128 ml m\ hk L[ - 32 ml m\ hk L^ 

-40 ml fXK + (16 ml + 160 ml m^) ^ik + (16 ml + 32 ml m\) fXK LI 




Bo{mu + md) = 2Bom, 



8{ml + 2m^)(2L^ - L^) + 8m^(2L^ - U,) + - i/i, (48) 



in agreement with |T6[. The 0{p^) contributions are 
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+ (-32 ml - 192 ml m\) fiK L; + (-32 m^ - 64 ml m\) fiK 

+ (-5/9 ml - 29/36 ml) /i^ + (32/3 m^ - 128/3 ml m\) jj^ L\ 

+ (8/3 ml - 32/3 ml m\) fir, + (8/3 ml - 32/3 ml m\) /i^ 

+ (-16 ml + 48 m^) /x^ + (-64/9 ml + 64/9 m^) //^ 

+ (64/3 - 160/3 ml m\) /i^ L^: + (-128/3 ml + 128/3 ml m\) fir, Ll^ 



-16/3 ml fir, LI + (-128 m^ - 512 - 512 m^ m\) 

+ (-256 ml - 384 m^ - 128 m^) LI LI 
+ (512 + 2048 ml m\ + 2048 ml m^) L;: 
+ (512 ml + 768 + 256 ml m^) 

-128 + (512 + 768 m\ + 256 m^) LI L\ 



r2 



+512 m° LI LI + (-512 m° - 2048 m^ - 2048 m^ m%) Lg 



r2 



.6 



1536 ml m\ - 512 ml m%) Uq LI - 512 ml Lf 



+ (-1024 m: 

+5/ 6 H^{ml, ml, ml; ml) ml -5/8 H^{ml, m^, m^; ml) ml 
+1/18 H^{ml, mj, mj; m^) + H^{m\, ml, m%; ml) ml m\ 
-5/6 H^{m\,m\,ml^]ml)ml - 1/8 i^^(mj, m^, m^; m^) 

I 1 /o rrF/ 2 2 2 2\ 2 2 , ttF f 2 2 2 2\ 4 

+ 1/2 H (m^, m^, m^; m^) m^ m^ + i^^ (m^, m^, m^; m^) 

I o jiri^/ 2 2 2 2\ 4 I o zjF / 2 2 2 2\ 4 

+2 (m^, m^, m^; m^) + 3 H^^im^, m^, m^; m^) 

o /o tji''/ 2 2 2 2\ 4 I o TjF / 2 2 2 2\ 4 

+9/8 i/fi(mj, m\,m]^]ml)ml. 

Where the definitions of H[ {m{, m^, m^; q^) appear in App. |C[ 
For the kaon we obtain 



(50) 



m. 



m\ 



K 



OK 



){4) 



B^ims + m) 



[m„ 



2m\){2Ll-Ll)+8ml{2Ll 



ml + 3mz 



6m^ 



-fir, 

(51) 



in agreement with |jT6[. The 0{p ) contributions are 



p2/ 2 \(6) 



- 16 m! Cl^ - 48 + 16 C\, + 48 C- 



K ^19 



+48 m\ C^^ + 48 m\ C^^ - 32 + 32 ^[^ 

-16 + 64 Cl^ - 32 C'l^ - 96 C[g 

-32 Cs'q + 192 ml + 32 ml Cg^s - 32 - 64 m\ Cl^ 

-32 Cl^ - 32 - 64 m% + 96 m% Cl^ + 128 m\ C^^ 

+192 Cs'i + 32 C?-, + 64 m% CI. 



(52) 



r4/ 2 \(6) 



K ) Loops 



(167r 



2^2 



1 /8 m\ vr^ - 13 /24 ml m\ - 73 /648 ml vr^ 
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-19/108 ml m\ - 763/1296 m\ tt^ - 4709/1728 m\ + -— - 3/2 m\\x^ 

+ (-3/2 ml m\ - 3/2 m]^) fiK + (-1/2 ml m\ - m%) fi^ + Am% L{ 
+ (56/9 ml m% - 16/9 ml m% + 122/9 m%) + (41/27 m^ m\ 

-A/21 ml + 89/27 m^) + (-2 m;' - 27/8 m^) /x^ - 3 //^ //^ 

+(1/3 ml - 41/12 m^) //^ //^ - 96 /x^ - 24 m\ /i^ 

—30 m\ //tt -^3 + (136 ml m\ + 64 m^) //^ L4 

+(24 + 32 m^) //^ + (-176 ml m\ - 128 m^) /i^ 

+(-48 ml m\ - 64 m^) + (-3/2 - 251/18 m^) /x^ - 8/3 m\ I^k I^v 

—144 m%iJ,K L\ — 72 m'^nx L2 — 60 mj^nx 

+(16 + 160 m%) hk L\ + 64 m]^ fix LI + (—32 ml m\ — 256 m^) jiK LI 

-128 + (-43/288 m^ m'^ - 25/32 - 5/9 m^) /tj 

+ (32/3 ml ml. - 128/3 m^) /x^ + (8/3 m\ - 32/3 m^) //^ L\ 

+ (14/3 - 56/3 m^) /x^ L^^ + (-8 ml m\ + 64 m^) //^ 

+ (16/9 - 8/3 ml m\ + 224/9 m^) /x^ + (16/3 ml m\ - 256/3 mi.) /x^ L^: 

+ (-64/3 + 64 ml m\ - 128/3 m^) /i^ + (-32/3 m^ + 112/3 m\ 

-64 m^) /i^ + (-128 m^ m\ - 512 ml m% - 512 m%) L^^ + (-128 m^ m\ 

-256 ml - 384 m^) LI + (512 m\ + 2048 rr4 + 2048 rr^) 

+ (128 ml ml + 384 ml + 1024 m^) LI L^ + (-64 m% - 64 rr4) 

+LI LI (256 + 512 ml m% + 768 m%) + (128 ml m% + 384 m%) 

+ (-512 - 2048 ml - 2048 m%) L^^ + (-256 

-768 - 2048 m%) LI L^ - 512 m% L^^ 

I o I A ZjF / 2 2 2 2\2 2|0/o ljF/ 2 2 2 2 \ 4 

+3/4 i/ {m^,m^,mj^;mj^) m^mj^ + S/8 H {m^,m^,mj^\m^)mj^ 

+ l/AH^{ml, m\, m\\m\) m]^ - 3/32 E^{m\, ml, ml; m%) m]^ 

+9/16i/^(m^,m^; m^; m^) + 3/4 H^{m\,m\,m\] m^) 

+181/288 if^(m^, m% mj; m^) - 3/2 H({ml, ml, m\; m\) 

-3/2 iff (m^,m^,m5;m^)m^ - 5/AH[{ml^,mli,ml;m%)m% 

+9/4 //fi(m2 , m^ mi.; m^J m^. - 9/32 i/fi(m^, ; m^ m^) 

+27/16 H^i{m\, ml, mj; m^) + 9/4 H^^{m\, m\, m\; m\) 

+27/32 H^Ml. K ml; ml) m% . (53) 

For the eta we obtain 
2 

"^or, = -Bo-(2ms + ?7i) , 

^.'K)^'^ = -8mJ(m^ + 2m^)L^ + ^mj(m^ - 4m^)L^5 
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16 16 



+ (-m^ + 2mi^)fiK + (^g< - ym^j /t^, + — » " (54) 
in agreement with |jT6|. The 0{p^) contributions are 



F^{ml)P^ = 32/27mlCl, - 32/9mlCl, + 16/3m^C[, - 16/9mlCl, 
+16 ml + 16/3 ml CI, + 128/9 ml Cl^ - 16 ml Cl^ + 16 C;^ 
-16 ^2^1 - 32/3 ml C;^ + 32 ml C^^ - 128/9 m^ + 64/3 m\ 
-320/9 m^^ Cl^ + 32/3 - 256/3 m^^ Cl^ - 320/9 CI, 

-256/3 mt C^g + 192 m^ m\ Cl^ + 64 m^ m\ C^^ + 128 m]^ C;^ 
-64/3 ml ml C;^ + 512/3 m^ m\ C;^ + 512/9 ml m^ + 640/9 ml m^ 
+320/3 ml m% + 640/9 ml m^ CI, + 128 ml m% ^[g - 384 ml m% C^g 
-192 ml ml C^q + 192 ml m^ C;^ - 256 ml m^ C^^ - 256/3 ml m^ C;^ 
-1024/3 m^m^C^3 - 2048/27m^C[2 - 1024/9 m^C[3 - 512/9 C^^ 
-512/9 ml C[5 - 256/3 ml C[g - 512/9 ml CI, - 512/9 ml ^[g + 256 ml Clg 
+256 C2"o + 256 ml C;^ + 512/3 ml C^^ + 512/3 ml C^^ + 512/3 Cg'g , 

(55) 

^'K)Sops = - 91/11664 m^vr^ - 1781/15552 - 269/486 vr^ 

-4133/1296 ml ml + 133/108 ml ml vr^ + 367/48 ml ml - 1091/729 ml vr^ 

-7567/972 m^) + ^ (^/i^ (2/3 - 4/3 ml ml - 16/3 ml) + (-4/27 

+ 16/9 ml ml - 64/9 ml ml + 256/27 m^) + (-58/27 + 88/9 ml ml 
-88/9 ml + 544/27 m^) + (-20/27 + 32/9 ml ml - 34/9 ml 

+ 152/27 m^)^ + /i^ (43/12 - 3m^) + /x^/i^^ (-8/3 - 32/3 m^) 

fx^n (-35/9 + 64/9 ml) + /i^ (32 ml - 128 ml) + /i^ (8 ml 
-32 m^) + jj^ L; (8 - 32 ml ml) + LI (-160/3 + 144 ml ml 
+256/3 m^) + fi^ LI (16/9 - 64/3 ml ml + 512/9 ml) 

Lq (64 — 544/3 ml ml — 512/3 ml) + /i^ (-384 ml + 1664/3 
-512/3 ml) + /i^ (-144 ml + 640/3 ml - 512/3 m^) 
+/i^ (-4 ml + 43/3 ^2 _ 212/9 m^^) + /i^ /i^ (28/9 ml - 64/9 m^) 
+fiK LI (128/3 - 512/3 ml) + L^ (32/3 - 128/3 ml) 

+fiK LI (56/3 ml ml - 224/3 m^) + fiK ^4 (-16/3 - 32/3 ml ml 
+256 ml) + /ix (80/9 ml - 32 + 1024/9 ml) 

+I2K LI (32/3 ml + 64/3 ml - 256 m^) + /i^ ^.'7 (-512/3 
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+1280/3 ml m% - 256 m%) + fxx LI (-224/3 + 704/3 ml m% - 256 m%) 

(-5/9 ml + 55/12 ml - 25/3 m^) + /x^ L[ (-16/3 + 128/3 ml ml 
-256/3 m^) + fXr, LI (-16/3 m^^ + 128/3 ml m\ - 256/3 m^) 
+//^ (-8/3 + 64/3 ml m\ - 128/3 m^) 
+At^ (32/3 ml - 112/3 + 128/3 m^) 

+//^ (m/27ml - 1472/27 m^m^ + 1408/27 m^) 
+fjLrj LI (-128/9 + 736/9 ml m\ - 2048/9 m^) 

(-1024/9 + 3584/9 ml m\ - 2560/9 m^) 

(-592/9 + 2176/9 rr4m| - 256 m^) 
+Lf (128/3 ml - 512 - 2048/3 m^) 

LI (256 /9 - 128/3 m^ - 128 - 5632 /9 m^) 

(-512/3 ml + 2048 rr4 + 8192/3 m^) 

(1024 - 3072 ml + 2048 m^) 
+LI LI (1024/3 - 256/3 m^ m\ - 3328/3 ml + 7168/3 m^) 
+Lf (-128/9 + 256/3 m^ m\ - 256/3 - 1024/9 m^) 

LI (-1024/9 ml + 256/3 m^ m\ + 1280/3 + 10240/9 m^) 

+L'5 (1024 - 4096/3 m^ m\ - 1024/3 + 2048/3 m^) 

LI (3584/9 - 6656/9 m^ m\ - 512/9 + 8192/9 m^) 

(512/3 ml - 2048 rr4 - 8192/3 m^) 

(-4096/3 + 4096 ml - 8192/3 m^) 
+LI LI (-1024/3 ml + 512/3 m^ m\ + 3584/3 ml m% - 4096 m^) 

(-4096/3 + 4096/3 m^ m\ + 4096/3 ml m% - 4096/3 m^) 
+Ll^ (-512 + 2048/3 m^ m\ + 2048/3 - 4096/3 m^) 

+1/6 H^{ml ml m^ m\) ml + 1/8 E^{ml, m\, m\, m^) ml 
+3/2 if^(m^, m\, m\\ mj) ml m\ - 1/2 E^{m\, m\, m\\ mj) 
+4 H^{ml^, m\, m^^; m^^) - 64/9 H^{m\, m^, m^^; m^^) 

I 1 /o tJ-F'/ 2 2 2 2\ 4 tltF/ 2 2 2 2\ 2 2 

+1/8 H {m^,mj^,mj^;m^)m^ - H [m^,m^,mj^]m^) m^mj^ 
+2 H^{m'^, m\, m\] m^) + 49/486 H^{m^^, m^^, m^; m^) ml 
-112/243 i/^(mj, mj, m^; mj) + 128/243 i/^(mj, mj, mj; m^) 

I TjFf 2 2 2 2\ 4 ^ zjFf 2 2 2 2\ 2 2 

+^1 ("^TT! "^X! "^x; "^r,) - 4 iii (m^, m;^, m^; m^^ 

I o LT-F/ 2 2 2 2\ 4 /o ZLT^'/ 2 2 2 2\ 2 2 

+2 (m^, m^, m^f, m,,) - 40/3 {m^, m^^, m^; m,^) m^ m^ 
+64/3 H[{m%, m\, m\\ mj) + 3/8 i?2^i(m^, rr^, m^; mj) 

o tJ-i^/' 2 2 2 2\ 2 2 I tjF / 2 2 2 2\ 4 

-3 i^ai ("^TT, "^x, "^i^ ; "^r,) + 6 H^iim^, m^, m^; m,,) 

I o /o ZjF / 2 2 2 2\ 4 o zltF / 2 2 2 2\ 2 2 

+3/8 H^i{m^, m^: 171^, m^) m^ - 3 H^i{m^, m^, 171^, m^) m^ m^ 
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+6H^^{ml,mj^,m%;ml)m%) . (56) 

A. 2 Decay Constants 

The decay constants are given by 

Fa = Fo{l + Fi'^ + (Fi^))cT + {Sloops) . (57) 

For the pion we obtain 

F^^Fi^) = A{ml + 2ml)Ll + 4m^L^ - 2/.. - /x^ (58) 
in agreement with |16[ for the 0{p'^) contribution. The 0{p^) contribution is 



Fl{F^f?T = 8 ml CI, + 8 CI, + 24 C[, + 8 C[, 

+16 m^, C[5 - 32 ml ml C{q + 32 , (59) 

Ft{Fl''^)ioops = (^35/288 + 41/128 + 1/144 m^m^vr^ 

-5/32m2 ml, + 11/72 m^vr^ + 15/327^1.) + -1- f/i, (1/2 + m^,) 



167r2 

+ (1/4 ml + m^) /iK + (1/3 ml - 1/3 m|) /i, - 2 + (-37/9 

+8/9 m^ - 52/9 m^) + (-28/27 m^ + 8/27 ml m\ - 43/27 m^) 
+ (8 + 20 + 8 m^) + (8 + 4 m^) + (-16 m^ 



^|+{^/% + m-^ml)^^l 



-40 - 16 m^) + (-16 - 8 m^J F 

+2/3 /i^ /i^ + 56 /ijr -^^1 + 32 L2 + 28 /i^ L3 

+48 /i^ LI + 12 /i^ + (-32 - 64 m\) /i^ - 32 ml 

+ (3 - 1/2 m^m^^) /i^ - 2/3 fiR fJ-q + 64 /ii^ + lQm\^K L'2 

+20 f^K LI + (12 - 8 m^) /ix ^4 + (-4 ml + 8 m\) ^ik LI 

+ (-16 - 32 m^) L^ - IQmj^ jjx L^ + 3/8/^^ 

+ (-16/3 ml + 64/3 m%) fi., L\ + (-4/3 ml + 16/3 m\) //^ 

+ (-4/3 + 16/3 m^) /i^ + (8/3 ml - 32/3 m\) jj^ LI - 4/3 ml /i^ Li; 

+ (-8 - 32 ml m], - 32 m^) + (-16 m^ - 32 m^) 

-8 + 5/12 H^\ml, m^ ml; ml) m\ - 1/2 {ml, ml, ml; ml) ml 

ri-ia ijFlf 2 2 2 2\ 4 I 1 /-la uF / 2 2 2 2\ 2 

-5/lQH {m^,mj^,mj^;mJm^ + l/lQH {m^,mj^,mj^;mj m^ 
+ 1/36 H^' {ml, ml^,ml; ml) ml + 1/2 H^'{mj^, ml, mj^; ml) ml ml^ 

1 /o ziri^/ 2 2 2 2\ 2 cr /i o zirF// 2 2 2 2\ 4 

-1/2 H {m^, m^,mK;m J mj^ - 5/12 H {mj^,mK,m^;mJ m^ 

1 /ic zir^''/ 2 2 2 2\ 4 I 1 tji^'/ 2 2 2 2\ 2 2 

-1/16 if (m^,m^,m;,;m^)m^ + 1/4 if {m^,mj^,mj^;mj m^mj^ 
+ 1/1Q {ml, m\,m\; ml) ml - 1/4 //^(mj, m^, m^; m^) 



25 



+1/2 HCiml ml, m^; ml) m% + El\m\, m\, mj; ml) mi 

I o /o TjFl/ 2 2 2 2\4 o l^r ZjFl/ 2 2 2 2\ 4 

+3/2 i/f/(m^, , m^; ml) + 9/16 //^/(mj, m^, m^; ) • (60) 

The functions H['{m\,m\,m\\q^) are defined in App. 0. 
For the kaon we obtain 

FlF^ = A{ml + 2ml)Ll + 4m^L^ - ^/i. - ^/i^ - \^^, • (61) 

The 0{p^) contribution is given by 

FI{fP)ct = 8 mt CI, + 24 - 8 m^ CI, - 16 ml ml CI, + %mlml CI, 

-32 ml ml C[g + 16 ml ml C\, + 16 ml CI, + 16 ml C[^ + 32 Cl^ , (62) 

i^.'(i^f )loops = (l/lQmU^ + 13/48 mt + l/32m>|,7r2 

-29/192 + 3/16 ml tt^ + 197/384 m]^ + (^15/16 /x^ 

+ (3/4 + 1/8 ml) fiK + (1/12 + 41/48 ml) /x^ - 2 mf^ L[ 
+ (-28/9 + 8/9 m^m^ - 61/9 m^) + (-41/54 m^ + 2/27m^ ml 
-89/54 m^) + (2 m^ + 14 ml ml + 20 ml) U, + (10/3 m^ - 8/3 ml ml 
+34/3 m^) LI + (-4 - 28 ml ml - 40 m^) + (-16 m^ + 32 ml ml 

-16 m]^) ^7 + (-12 mt + 16 ml ml - 28 ml) L^^ + (9/32 + 3/8 m'^ ml) 

+3/8/1,^/1 A- + (13/16 + 3/4 m^ m;j^^)/i^/i^ + 48m^/x^ L\ + 12 m^ /i^r ^2 

+15 ml /i^ L3 + (-15 ml + 18 ml-) yU^ L4 + (6 m^ - 3 ml-) yU^ L'5 

+ (-12 ml - 24 m^) /x, - 12 m^ + (27/8 + 3/4 ml m^^) 

+ (-29/8 - 1/2 m^m~^) /ix /"r? + 72ml fix + 3Qml fix 

+30 m^ /ii^ L3 + (18 m^ + 4 ml) fix Ll + Q ml fix LI 

+ (-24 ml - 48 m^) fix ^ - 24 m^ fix ^ + (17/8 + 1/32 ml m'^) /ij 

+ (-16/3 ml + 64/3 m^) /i^ L\ + (-4/3 m^ + 16/3 ml) fi^ 

+ (-7/3 ml + 28/3 m^) /x^ + (35/3 m^ + 22/3 ml) /i^ 

+ (-2 ml + 23/3 m^) LI + (-12 m^ - 24 ml) /i^ + (-18 m^ m-^ 

+42 m^ - 24 ml) /x^ + (-6 m^ m'^ + 18 m^ - 24 ml) /x^ + (-8 m^ 

-32 m^ m^ - 32 ml) Lf + (-16 ml ml - 32 ml) LI LI -8 ml Lf 

+3/8 H^'{ml, ml, ml; ml) ml ml + 3/16 H^'{ml, ml, ml; ml) ml 

-3/8 H^'iml, ml, ml; ml) ml + 1/8 H^'{ml, ml, mj; ml) ml 

-3/64 H^'iml, ml, ml; ml) ml + 3/64 H^iml, ml, ml; ml) ml 
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+9/32 H^'{ml, ml, mj; ml) m% - 9/32 H^{ml, ml, m\- m\) m\ 

I o /o TjFn 2 2 2 2\4 o/o ZLT^'/ 2 2 2 2 \ 2 

+3/8 -H {mj^,mj^,mj^;mj^j mj^ — S/8 H [mj^,mj^,m^;mj^) m^^ 

+181/576 H^\m\, m^, m^; m\) m% - 9/64 H^{ml, m^, m^; m\) m\ 

~?,/AH('{ml,ml,m\;m\)m]^ - ?>/AH(\m\,ml,ml;m\)m% 

-5/8 iff (m^, mj, mj; m|-) + 9/8 i^af "^x ) "^/c 

-9/64 iif/ (m^, , ml; m^) m% + 27/32 iif (m^, , m^; m^) 

+9/8 H^^iml, m\, m\; m\) m% + 27/64 H^^{ml, m^, m^; m\) m% . (63) 

For the eta we obtain 

F^^Ff = A{ml + 2ml)Ll + ^(4mi - m^)^^ - 3/.^ (64) 



in agreement with |T6[ for the 0(p ) contribution. The 0{p ) contribution is 

F^{F^'^)cT = 8 mt CI, - 8/3 mt CI, + 24 C[e + 8 + 64/3 

-64/3 ml m\ CI, + 16/3 ml m\ CI, - 32 ml m\ - 64/3 ml m\ C{j 

-128/3 ml m\ C^g + 64/3 + 64/3 CI, + 32 C[g 

+64/3 m], CI, + 64/3 , (65) 



F,^(F^(^))loops = (^5/96 mtvr^ + 91/384 + l/48m2 m^vr^ 

-1 1/96 m^^ + 5/24 m| vr^ + 49/96 m^) + ^ (^(3/4 + 3 m^) /Xi^ 

-m|/i^ + (-2/9 mt + 16/9 32/9 m^K)L[ + (-29/9 m^ 

+ 16/9 - 68/9 m%) + (-10/9 m^ + 8/9 - 19/9 m^) 

+ (12 ml m\ + 24 m^) + 12 + (-24 ml m\ - 48 m^) 

-24 + 9/8 /x^ + 48 ml fi^ L\ + 12 ml /x^ + 12 ml /x^ 

-24 /i^ - 4m^/i^L^ + (3 + 3/2 rn]f) fxj^ - 2/ix/ir, 

+64 /iii- + 16 mj^ hk L'2 + 28 mj^ hk + (36 ml + 40 m^-) yU/^ L4 

+ (4 + 56/3 ml) fiK L; + (-48 ml - 96 m^) /i;^ L^-ASmluK 

+9/8 /ij + (-8 ml + 32 m^^) /i^ + (-8 + 32 m^) fx^ 

+ (-4 ml + 16 m^) /i^ L^^ + (8/3 - 32/3 m^) /i^ + (28/9 ml 

-64/9 m^^) fi^ LI + (-8 - 32 ml m\ - 32 m^) + (16/3 m^ 

-32/3 - 128/3 m^) + (-8/9 m^ + 64/9 m\ - 128/9 m^) 

+ 1/12 H^\ml, ml, mj; m^) + 1/16 //^'(m^ , m^, m^; mj) 

\ o I A zjFi f 2 2 2 2\ 2 2 n/irzirf'/ 2 2 2 2\ 2 

+3/4 if [m^,mK,mK;m^)m^mj^ {m^,mj^,mj^;m^) m^ 

1 //I O'i'''^ 2 2 2 2\ 4 I r, TirF// 2 2 2 2\ 2 2 

-1/4 ii {mj^,mj^,m ■,m )m^ + 2H {mj^,mj^,m ;m )m^mj^ 
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— 32/9 i?^'(m|-, m|^, m^; m^) + 1/16 H^^ni^, m\, mj^; m^) m\ 

1 /o U-F'/' 2 2 2 2\ 2 2 I rrF// 2 2 2 2\ 4 

-1/2// {m^,mK:mK;m^)m^mK + H {171^,171^,171^,171^) 171^ 
+3/16 //^(mj, m^, m^; mj) - 3/4 //^(mj, ml, ml; m^) 
+49/972 H^\ml, mj, mj; mj) - 56/243 H^\ml, mj, mj; m^) ml ml 
+64/243 //^'(m?, K ^S) + 1/2 i^f V' , ^x, ^x! ^J) < 

o tJ-i^'/' 2 2 2 2\ 2 2 I TjFif 2 2 2 2\ 4 

-2-f^i {m^,mK,mK]m^)m^mK + H-y {mj^,mj^,m^;m^) m^ 
-20/3 H['{ml, ml, m\; mj) + 32/3 H['{ml, ml, m\; mj) 

+3/16 H^l{ml, ml, ml; mj) - 3/2 H[-^{ml, ml, ml; mj) ml 

I o Eji^'/ 2 2 2 2\ 4 I o /ir rj'f''/ 2 2 2 2\ 4 

+3 ii2i ("^TT, "^x; "^,?) "^x + 3/16 H^i (m^, m^^, m^^; m^) m^ 

-3/2 i/f/(m;, m^, ml; mj) + 3 i/f/(m;, m^, ml; mj) . (66) 

B Axial- Vector Two-Point Functions 

B.l The remainder transverse part 

The remaining transverse pieces are given by 



CT \ ^U'^-^ W loops 



a=Tr,K,r] 



(67) 



The O(p^) contribution is the same for all three two-point functions 

{m(s')f = {m,'))'" = = 4L;„ -SHI. (68) 

The (!?(p®) contribution can be written in a somewhat simpler fashion by using 
the function 

TjT / 2 2 2 2 2\ o AttFi/ 2 2 2 2\ 2 ttF / 2 2 2 2\ /cnN 

ii [mi,m^,m^;q ;m^) = Zm^H^^[m-^,m2,m^;m^) - q H^i[mi,m^,m^;q ) .(69) 
The isospin remainder is 

,(6) 



(g2fl3i(g2))^; = -8 C^^, - 32 ml C{, ^ - 32 ml C^^ ^ - 32 ml C;, q' 
-32 ml CI, - 64 ml g^ - 64 ml CI, g^ - 16 CI, q' - 8 C^, q' , (70) 



- (I6^ (3/32< + l/6<.^.^ + l/3m^g^ 
+1/12 ml q^ TT^ - 5/24 ml q^ + 3/32 g^) + (2 ?'/^^ + ) 
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+5/3mlH^\ml,ml,ml-,ml) - 2 ml {ml, ml, ml; ml) 
+2 ml H^{ml, m^, m^; q') - 5/4 m^ H^'{ml, m\, m^; ml) 
+1/4 ml H^{ml, m\, m\; ml) - 1/4 ml H^{ml, m\, m\; q^) 

I 1 /n 4 TjFlf 2 2 2 2\iri2 2 rrF// 2 2 2 2\ 

+1/9 m^// {m^,m^,m^;mj +2m^m;^// (mK,m^,mK]m^) 

o 2 tj-F/ 2 2 2 2n I o 2 rrF/ 2 2 2 2\ 

-2mKH (m^,m^,m^;m^) +2m^i/ {mK,m^,mK;q ) 

r- /o 4 tltF// 2 2 2 2\ I / 2 tjf''/ 2 2 2 2\ 

-5/3 m^ii [m^,mj^,m^;m^) + [m^H {m^,mK,mj^;mJ 

-H^iml, ml, ml; ml) + //^(mj, m^, m^; g^)) (-1/4 ml + m^) 

+2 ml H['(ml, ml, ml; ml) + 4 H['(ml, ml, mj; m^) 
+2 //^(m^, m^, m^; g'; ml) - 1/4 //^((m^, m^, ml; q^; ml) 
+2 //^(m^, ml, ml; q^; ml) + 3/4 //^(mj, m^,, m^; g"; ml) . 

The remainder for the kaon two-point function is 

(?'nW = -32 CI, q' - 32 C;, q'-Sml C;, - 32 ml CI, 

-64 C[3 g^ - 32 - 64 g^ - 16 g^ - 8 C;, q' , 

(fFlii^AK{<f))Z,, = JY^, (l/24 ml q' ^ - 5/48 ml q' + 3/32 ml 
+5/24 ml g2 tt^ + 11/48 ml q^ + 3/32 g^) + ^ ^3/4 g^. 

+3/2 q^iiK + 3/4 g Vr, ) + 3/2 m;' g^' - 12 gV- 

+3 ^2^2^ _ 24 2,?;^ ^2^^ ^ ^-2 ^2^2 _ ^2^^ 

+3/ 4 ml (2 + m^) H^\ml, ml, ml; ml) 

0/0 2 ttF/ 2 2 2 2 \ 

-3/2 m^H {m^,m^,mK;mK) 

+3/2 ff^K , , m^; g^) + 1/2 i/^'K , ml, m\; ml) 
-3/16 ml H^'{ml, ml, ml; ml) + 3/16 ml H^{ml, ml, ml; ml) 
-3/16 ml H^iml, ml, ml; q") + 9/8 ml H^'{ml, ml, mj; ml) 
-9/8 ml H^{ml, ml, mj; ml) + 9/8 ml H^{ml, ml, mj; q^) 
+3/2 H^'{ml, ml, ml; ml) — 3/2 m|- H^{ml, ml, ml; ml) 
+3 / 2 ml {ml, ml, ml; q^) + 181/144 /J^'(m^, mj, mj; m^) 
— 9/16m|- H^(ml, m^^, m^; ml) + 9/16 m|- H^{ml, m^^, m^; g^) 

o 4 rrF// 2 2 2 2 \ o 4 rrF// 2 2 2 2 \ 

-3mj^H^ {m^,m^,mj^;mj^) -Smj^H^ {mj^,m^,m^;mj^) 

r- /o 4 ttFu 2 2 2 2 \ I o /<-) ttT / 2 2 2 2 2 \ 

-5/2 m^iii [mj^,m^,m.^f,mj^) +3/2H {m^,m^,mj^;q ; m^^) 
-3/16 i/^(m^, m^, m^; g^; ml) + 9/8 i/^(rr4, , mj; g^; m^) 
+3/2 H^{ml, ml, ml; q^; ml) + 9/16 H^{ml, mj, mj; g^; ml) . 
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And finally the remainder for the hypercharge is 

,(6) 



+64/3 ml m\ C;, + 32/3 ml C^^ - 32 ml C^^ + 32/3 ml C;^ q^ 

-32 ml C;^ g2 ^ 256/3 C^^ - 64/3 C^^ - 128/3 q^ 

-64 C[3 g2 - 128/3 mj, C;, q' - 64 mj, C^, q' - 16 C^, q' - 8 , 

(74) 

i^"F^<i^"))Zs = (T^ (V96rnt - 1/12 - 1/4 m^g^ 

+ l/6rr4 + 1/4 m^g^TT^ + 3/8 rr^g^^ + 3/32 g^) + -^3gVi^ 



167r2 

+6 m^^ q^ - 48 L^q q^ /ix + 1/3 H^'{ml, ml, m^; m^) 
+ K + 9/4 m;) iJ^'K, m^, m^; ml,) 
-9/4 i/^(m2 , ml, ml; ml,) + 9/4 ml H^(ml, ml, ml; q^) 
+(-m^ + 8 - 128/9 + 9/4 mj) H^'{ml, ml, m\; mj) 

-9/4 {H''{m%,ml,ml;ml)-H^{m%,ml,ml;q^)^ 

+ (49/243 - 224/243 ml ml + 256/243 ml) H^'im^, m\, mj; mj) 

c 2 2 rrF// 2 2 2 2\ 

-Gm^m^H^ {m^,mK,mK;m^) 

\ A ( o 2,0 2\ 2 TjFl i 2 2 2 2\ 

+4(-3m^ + 8 m^jm^i^i [mj^,mj^,m^;m^^ 

+9/4 i/^(m2 , m^, g"; mj) + 9/4 {m\, ml, ml; q^; m\) . (75) 

B.2 The remainder longitudinal part 

In that case the 0{p^) contribution vanishes so we have that 

ftStf)^Ktf)):::+Kitf));:j_^,. (76) 

The expansion of the resummed selfenergy around the relevant pseudoscalar 
mass leads in general to rather high derivatives and produces naturally the com- 
binations 

ttLi 2 2 2 2 2\ 1 ( 2ttFi 2 2 2 2\ 

Hi {mi,m2,m^;q ;m^) = _ ^2^2 (^^ (^i'^2>^3;g ) 

2ttF/ 2 2 2 2\ 2/2 2\ ttF/ / 2 2 2 2\\ /^^\ 

-q {m^,m2,m^;m^) -m^{q - m^)Hi {m^,m2,m^;mi) j (77) 

ttM f 2 2 2 2 2\ 1 ( ttF I 2 2 2 2\ 

R {m^,m2,m^;q ;m^) = --^ IWl"^ ("^i, "^2> "^3; 5 ) 
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—H {rrii, m^; m^) — {q — ni^)H '{rn-^,m2,m^;m4)j. (78 
All of these functions are regular at q"^ = 



4- 



The longitudinal isospin remainder is 

,(6) 



(M°i(a)^; = 8m^C,^, (79) 
(6) -3/32 



loops (IGtT^)^ 

I 1 /o ZjL/ 2 2 2 2 2\ 4 o uM / 2 2 2 2 2\ 6 
I jltL/ 2 2 2 2 2\ 4 I 1 //I ttM / 2 2 2 2 2\ 6 

+i7 (m^, m^, m^; g ; m^) + 1/4 (m^, m^, m^; g ; m^) 

1 /n ziri/ 2 2 2 2 2\ 4 o ttM / 2 2 2 2 2\ 4 2 

-l/9if {m^,m^,m^]q ]m^)m^-2H [mK^m^^m^^q ]m^) m^nij^ 

I cr /o ziri/ 2 2 2 2 2\ 4 o //( zirM/ 2 2 2 2 2\ 4 2 

+5/3H (m^,m^,m^;g ;m^)m^ - 3/4 ii [m^.m^.m^^q ]m^)m^m^ 

o rr-L/ 2 2 2 2 2 \ 4 /i rrL/ 2 2 2 2 2\ 4 

-2 (m^, m^, m^; g ; m^) - 4 //^ (m^, m^, m^; g ; 

r ttL / 2 2 2 2 2\ 4io//(rrL/'2 2 2 2 2\ 4 

-QH2i{m\,ml,m\]q^]ml)ml - 9/4 i/2^i(mj, m|, m^; g^; m^) . (80) 
The kaon remainder is 

(5^n?ka)^'i = 8m^^C,^, (81) 



(6) -3/32 



loops (IGvr^)'^ 

o /"-) zirA^/ 2 2 2 2 2\ 2 4,0 /ic rjil^/ 2 2 2 2 2 \ 6 

-3/2 H (m^, m^, m^; g ; m^^) + 3/16 (m^^, m^, m^; q ; m^) 

-3/4 H^{ml, ml, m\] g^; rr4) rr4 - 1/2 H^{ml, m\, mj; g^; m^) m% 
-9/8 H^{m%, ml, mj; g^; m^^) - 3/2 H'^{m\, m\, m\; q^; m\) m% 
-25/36 H^{ml, mj, mf^; g^; m^) - 9/16 H^{m\, mj, mf^; g^; m^) 

I o Jir-L/' 2 2 2 2 2\4io zji/ 2 2 2 2 2 \ 4 

+3 [m^, m^, m^, q ; m^) m^ + 3 H-^ (jn^, m^, m^; q ; m^) m^ 

I cr /o zjLf 2 2 22 2\4 n/o ziri /2 2 2 2 2\4 

+5/2 (mj^, m^, m^; g ; m^) - 9/2H2i{m^, m^, m^, q ; m^) 
+9/16 H^i{m\, ml, ml; q^; m\) - 27/8H^^{m%, ml, mj; g^; m^) 

n /o zltL / 2 2 2 2 2 \ 4 ot/ic iiri /' 2 2 2 2 2 \ 4 

-9/2 H2i{mK,mK,mK;q ■,mK) - 27/16 H2i{mK,m^,m^; q ■,mK)mK. 

(82) 

Finally the hypercharge remainder is 

(g2<(g2))^'^ = -256/3 mtC3^, + 8 m^Cg^ + 512/3 m>|C3^, 

-64/3 ml m\ C^^ - 256/3 C^^ + 64/3 C^^ , (83) 
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1 /o ttL/ 2 2 2 2 2\ 4 ttL/ 2 2 2 2 2\ 4 

-1/3 H {m^,m^,m^;q ;m^)m^-H {m^,mj^,mj^; q ]m^)m^ 

n I A ttM/ 2 2 2 2 2\24 

-9/4 ii (m^, m^, m^; g ; m^j m^m^ 

+H^{ml, ml, ml; q^; mj) (m^ - 8 + 128/9 m^) 

-9/4/^^(771^, rr4, rr4; g^; m^) mf^ 

+1/243 //^(mj, mj, mj; q"^; mJ)(-49 + 224 ml m\ - 256 m^) 

+6i/f (m^, m^, m^; g^; mjm^mj + H^(mj^, m^, mj; g^; m54(3m^ - 8m^)mJ 

-27/4 //2'l("^' , - 27/4 //2^i(mJ, m^, m^; g^; mj)mj . (84) 

C Loop integrals 

We use dimensional regularization here throughout in d dimensions with d — 
4-2e. 

C.l One- loop integrals 

We need integrals with one, two and three propagators in principle. These we 
define by 

d^g 1 
(27ryq 

We also use below 



A{m') = 1 / ^ ^ ^ ■ (85) 

^ ^ i y 27rY q^-m^ ^ ' 



A{n, m^) = - f , ^ ^ , (86) 



i J (27r)^ (g2 - m2) 

which can be obtained by derivation w.r.t. of Aim?'). 
The two propagator integrals we encounter are 



B{ml,ml,p^) = T y" 
Bi,{ml,ml,p) -- 



d'^q 1 



{2T:Y{q^-ml){{q-pf-ml) ' 
1 /■ d'^g g^ 



B^^{ml,ml,p) ^ -J 



{2ttY (g2 — mi)((g — p^ — ml) 

D / 2 2 2\ 

Pf,Biim-^,m2,p ) , 
1 f d'^q q^^qi, 



{27iY (g2 — mf)((g — p)^ — m^) 
= Pij,PvB2i{ml,ml,p^) + g^yB22{rnl,ml,p^) . (87) 

All the cases with three propagator integrals that show up can be absorbed 
into the two-propagator ones by moving to the real masses rather than the lowest 
order masses. This provided in fact a consistency check on the calculations. 
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The explicit expressions are well known 

-Cln(m2))| + 0(e2), 

/ 9 9 9, 1 /, m? ln(m?) — mn ln(mn)\ -,9 9 9, , 

167r^ Y mf — m2 / 

J{m^,m2,p) = dx\n\ ^ — ■ ^ , 

levr"' JO \ mfx + m2{l — X) J 

C = ln(47r) + 1 — 7 and Aq = 1/e + C. The function J{rn\,m\,p^) develops an 
imaginary part for > (mi + 7712)^. Using A = — m^, S = ml + and 
1/2 — p4 _|_ ^4 _|_ ^4 _ 2p'^m\ — 2p'^m\ — 2m\m\ it is given by 

The two-propagator integrals can all be reduced to B and A via 
Bi{ml,ml,p^) = -^(^A{ml)-A{ml) + {ml-ml-p'^)B{ml,ml,p^)^, 

B22{ml,ml,p^) = ^ 1) (^("^2) + 2mlB{ml, mj, p^) 

i/2 2 2\ D / 2 2 2\\ 

+ (m2 - -p )Bi{m^,m2,p ) ) , 
52i(m?,m^,/) = ^(^A{ml) + mlB{ml,mlp^)-dB22{ml,mlp^)) . (90) 



The basic method used here is the one from Passarino and Veltman 25 



C.2 Sunset Integrals 

In this appendix we discuss the nontrivial two-loop integrals that show up in this 
calculation. They have been treated in several places already, in general and for 
various special cases. We use here a method that is a hybrid of various other 
approaches. We only cite the literature actually used. We define 

_ 1 r d'^q d'^r X 

J {2txY {2'kY (g2 - m\) {r^ - ml) {{q + r- pf - mf) ' 

H{mlmlml;p^) = ((!)), 

H^{m\,ml,ml]p^) = {{q^)) = p^Hi{ml,ml,ml;p'^) , 
H,Mlmlmiy) = {{q,q^)) 

TT / 2 2 2 2\| ZJ f 2 2 2 2\ 

= p^pyH2i[m^,m2,m^]p ) + g^^H22{m^,m2,m^;p ). 

(92) 
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By redefining momenta the others can be simply related to the above three. In 
particular 

{{%r^)) = Pf,p^H23{ml,ml,ml-,p'^)+ g^^H2i{m\,ml,ml-,p'^) , (93) 

with 

2H2i{ml,ml,ml]p'^) = -H2i{ml,ml,ml;p^) - H2i{rril,ml,ml;p'^) 

, TT f 2 2 2 2\ , nzj f 2 2 2 2\ 

+H2i[m^,m^,m^;p ) + 2Hi[m^,m^,m^]p ) 
+2Hi{rn\, m\, m\; p^) — H{rn\, m\, m\; p^) , 

nzj ( 2 2 2 2\ ZJ f 2 2 2 2\ tj I 2 2 2 2\ 

2E2^[m^,m^,m^\p ) = -H22{m^,m2,m^;p ) - H22{rn^,rni,m^\p ) 

+H22{mlml,ml;p^). (94) 

The first two follow from interchanging q and r and the third from the fact that 
it is proportional to g/^i, or p^Pv, which are both symmetric in /i and u. The last 
one follows from 

{qnTu + r^^Qu) = {Qiu + - Piu) {Qu + ru- Pu) - q,j,qu - r^r^ - p^p^ 

+2p^ {q, + r,) + 2p, {q^ + r^) (95) 

and redefining momenta and masses on the r.h.s.. In addition we have the relation 

p^H2i{m\,m\,m\-,p'^) + dH22{'i^\,'i^\,'rn\',p^) — 

mlH{ml, ml, ml\p^) + A{ml)A{ml) . (96) 

which allows to express H22 in a simple way in terms of H21. There is also a 
relation between and H 

Hi{m\,ml,ml]p^) + Hi{ml,ml,ml-,p^) + Hi{ml,mj,ml;p'^) = 

H {ml, ml, ml; p^), (97) 

which allows to write Hi{m^,m'^,m'^;p^) = 1/3 H {m^ , m^ , m^ ; p^) in the case of 
equal masses. The function H is fully symmetric in ml, ml and m|, while Hi, 
H21 and H22 are symmetric under the interchange of m| and m|. 

We do not exphcitly evaluate the integrals analytically. H, Hi and H21 are all 
finite after two subtractions. We therefore evaluate them as follows {Hi stands 
for H, Hi and H21) 

Hi{ml,ml,ml;p'^) = Hi{ml,ml,ml;0) + p'^-^Hi{ml,ml,ml;0) 

+Hi{ml,ml,ml;p^) . (98) 
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The functions Hi are finite in 4 dimensions and can be evaluated by their disper- 
sive representation |^6|, ^ or below threshold by the methods of . 

The value at zero and its derivative there have been derived essentially using 
the methods of except that we use a slightly simpler procedure than the 
recursion relations given there. 

First we define the intermediate integrals 



/(rai,n2,n3) 



2 J {2t:Y {2txY (g2 _ m\f^ {r^ - mlY^ {{q + - m. 



(99) 



which show up in the momentum expansion of Hi. The I{ni,n2,n^) with one of 
the rii = are separable and are e.g. 



/(ni, n2, 0) = A{ni, ■ml)A{n2, ml) . 



(100) 



All the others can be derived by taking derivatives of /(1, 1, 1) w.r.t. the masses 
ml, m^ and m^. The function /(1, 1, 1) is taken from p8|], note that our definition 
differs by overall factors from theirs. 



7(1,1,1 



-1 



r^fi 



-{An 



,2e , 



(167r2)2 2(l-e)(l-2e) 

2 ("^1 ~'~ "^2 + "^3) "I {^1 I'^l +"^2 1^2 +^^3 I1I3 

+ml (In2ln3-lni(lni + ln2 + ln3)) (ln3 lui - ln2(lni + ln2 + ln3)) 
+ml (Ini ln2 - ln3(lni + ln2 + ln3)) + ^(m^, ml, ml) + 0{e) \ . (101) 



In ( |101| ) we used luj = ln(m2) and the function '${ml,ml,ml). The expression 



for \E' is somewhat dependent on the relation between the various masses. Using 
Am = A(m^, ml, ml) = {ml — ml — mg)^ — Amlml , (102) 



we have for the case < 28 



^(m^^, m^, m^) = 2^— {Cl2(2 arccos zi) + Cl2(2 arccos2;2) + Cl2(2 arccos2;3)} , 

(103) 



with 

Zl 



—m\ + ml + m| 



Z2 



—ml + m| + m\ 



2m2m^ 2m^mi 
The case mi + m2 < m^, with Am > 0, is 



^3 



—ml + m\ + ml 



2mim2 



(104) 



\E'(m^, ml, m 



, m? ml 

' Am i 2 Inxi lna;2 — In — 5- m — ^ 



mn 



mn 



TV 



+ y -2Li2(Xl) 



2Li2(s2) 



(105) 
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with 



'-"a ' '-"1 '-"2 V '^m _ •■-a ' •■"2 V I'lnfi'l 
^1 - 7^2 ' ^2 - • (iUbj 



The cases mi + m3 < m2 and m2 + < nii can be obtained from the last one 
by relabelhng masses. Li2(a;) is the dilogarithm defined by 

U2{x) = - / -^ln(l-xt), (107) 
and Cl2(x) is Clausen's function defined by 

Ch{x) = -j^ dtln|2sin-| = -z(Li2(e") -Li2(l)) + (^O^) 

Notice that \E'(mi, m^, ^.3) is fully symmetric w.r.t. the three masses. The 
I{ni,n2, ris) for general rii can be obtained by taking derivatives of /(1, 1, 1). The 
relation 



(109) 

allows an easy evaluation of all needed derivatives and is equivalent to the recur- 
sion relations used in P^ . 

In order to express the Hi functions at zero and the derivatives w.r.t. 
at zero the easiest is to shift momenta to g = g — p in the integral and then 
Taylor-expand using 

I = -J— y (zltJLlZ]' . (110) 

{q + pY - ml g2 _ ^2 y - ml J 

The integrals can then be done using q^q^, — > Qtiu^'^/d and equivalent identities 
for the higher orders. We have run this procedure to higher orders then necessary 
to check the cancellations of infinities there. This results in 

H{mlmlml;0) = 1(1,1,1), 

■^H{ml,ml,ml-0) = i^/(2, 1, 1) + ^m?/(3, 1, 1) , 

H,{mlmlml;0) = ^J(l, 1, 1) - ^771^/(2, 1, 1)) , 
^H,{mlmlml;0) = ^^^i_^((2 - rf)(rf - 4)/(2, 1, 1) 

+8{d - 4)m2/(3, 1, 1) - 24mt/(4, 1, i; 
1 



H2,{mlmlml;0) = ^^^_-^(t/(d - 2)/(l, 1, 1) - 4(rf - 2)m2/(2, 1, 1) 
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+8mt7(3,l,l)), 

-^H2i{mlmlml-Q) = + lod _ 48)7(2, 1, 1) 

192 / 

+12(d - 10)m?/(3, 1, 1) - 72mt/(4, 1,1) + ^ (/(2, 1, 1) 
+3m^/(3, 1,1) + 3mt/(4, 1, 1) + m^/(5, 1,1))}. (HI) 
Evaluating these expressions then leads to 

{lQT:''fH{ml,mlml-Q) = ^A2(m? + ml + ml) + ^Ai(m?(l - 21ni) 

+ml{l - 2 Ins) + ml{l - 2 Ing)) - ^^(m?, m^, m^) + + ^ _ in, 

+ ^ ( - ln2 ln3 + Ini ln4 )^ + ml + ^ - ln2 +^ - Ini Ins + ln4 )^ 

+ml + ^ - Ing +^ ( - Ini ln2 + Ing ln4 )^ + 0{e) , (112) 

(167r^)^^i/(m?, m^, m^; 0) = -^A^ + ^^^H!^*(m?, m^, ml) + 1 

Im +7^(m^ -m\- ml) In^ +^{ml -ml- ml) Ing +C»(e) , 

(113) 



(167r2)2i/,(m2,m2,m2; 0) = 'I^^±^\^ + lAi(2m2 + ml{l - A\n^) 

4 8 ^ 

+m^(l - 4 Ins)) + ^ ( - 1 + ^) *(m?, ml ml) + - ^ 



2/^ 9 1, 1/, n 

+^2 I ^ + ~ i +^(ln2 ln4 - Ini Ing; 
9 1, 1 



+""'1 24 + 16 " i " ^""'^ ' + ^^'^ ' ^^^^^ 



nfi^2^2 u (^2 ^2 ^2.n^ _ ^ \ , ^1^2^3'^m ,T,^^2 ^2 ^2\ ^l^rn 

^ i Q^^^y^i^^2: "^3> Uj - -^2^ A3 W(mi, mgj - 

mflni/ mamo\ ln2 / 90 9 6m?m? _ 9,\ 

+^ (^3 - - 2m? - ^^(5^ + 2ml)^ + ^ + ^(e) , (115) 
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(167r")"if2iK,m^,m^;0) 



m2 + rrig 



A2 + T^Ai 1 3ml + ml{2- 12 Ins) 
6b \ 



+m2(2- 121n3) + 



+ml\ 



17 1 

liii 4- 

72 3 



6A,l V 

n 



6A, 



19 



1 In 



m 



+ml 



3Q 54 9 

19 1 

\ Ing - 

36 54 9 6A 



-{6m + 2m\ ) ln2 + - (Ins ln4 - Ini In, 



TD^ 1 \ 

^(5m + 2m^)ln3+-(ln3ln4-lniln2) + C>(e) , 

(116) 



(167r2)2 AiJ2i(m2 m^, m^; 0) = -— Ai + — 
^ ' dp^ ^ ^ ^ ^ ' 24 288 



m 



24A, 



'-{6m + 2ml) 



6Am 



fi 2 2 

m{m2m^ 



Am 



mimr,m^ 



1 



In^ 



dm'^m'^ 



A. 



\l/(m^, mg) + 



I'-'m 



12A^ 



Ini 1 + 30- 



m%m% 



ml — ml) + 



m% 



m 
.2^2 



— — In3(m3 - ma) 

iZAm, 



\t 2 ■ 

ln3 — 3m;L — 12mim2 + 5m]^m3 + 3m2 — ^2^3 — 2m3 



12A^ 

4 2 4 

m\m^m2 



A3 

*m 

2 2 

m|m2 



ln2-(m; 



12A2 



Ino 



Here we used 



mo 



+ 



ln2 (m: 



12Am 

3m^ — 12mlm\ + hmlm\ + 3m3 



m.. 



m^mn 



2m% 



0{e). 
(117) 



ln4 


= Ini + ln2 + ln3 , 




6m 


222 
= m^ — m2 — mg , 




A2 


= Ag + (ln(47r) + 1 




Ai 


= Ao + ln(47r) + 1 - 


-7, 


Ao 


= - + ln(47r) + 1 - 


7- 



(118) 



These are appropriate for MS subtraction as is customary in CHPT 29 



Below threshold the methods of p^l l^ad to a two-integral representation of 
the finite part 

{H, Hi, H2i}{ml, ml, m\; f) = 
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/ rf(T. A 1,— ,— X / (ix/C2(x,a,p^){l,x,x^} , (119) 

J(m2+m3)2 \ \ a a J Jo 

with 

_ . 2\ 1 /i ?^?(1 — a;) + cx — xfl — p^xll — x) \ 

lC2{x,a,p ) = ——^ \\n J- — + 



(IGtt^)^ y mf (1 — x) + ax mf (1 — x) + ax J 

(120) 

and 

A(a:,?/,2;) = (x-y-2;)^ -4?/2;, (121) 

the Kallen function. 

The dispersive representation 



d-3 



B{mimit) = —-j—^ i rfa A 1,-^,-^ a^ 



2 



5 



(122) 

has been used here instead of the simpler case with equal masses used in p7[] . 

Above threshold, s > (mi + m2 + m^Y functions develop imaginary 
parts and they can then be evaluated from their dispersive representation 



Hi[m^,m^,m^;s) = — — . (123) 

TT J (m,i+m2+m3)2 Z Z — S 

The imaginary parts are given by (in d = 4) 



Im Hi, if2i} {ml, ml, ml; s) = 



with 



= -^{s + ml-im^ + m^f), 

l\J s 



^23 



Tpmax jpmin 
^2 ~ ^2 



s + ml- 2^Ei , 

2 r ^Ks,mi,ml^)^\{ml^,ml,ml). (125) 

77123 V"S 



In the text we use mainly the finite functions H[ {m\,m\,m\\q^) and their 
derivatives with respect to g^, H['{m\,m\,m\\q^). These correspond to using 
Eq. (^) and setting A2 and Ai to zero in Eqs. (|112| - |117| ). 
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D Regularization and Renormalization 



In this paper we have employed the version of Modified Minimal Subtraction 
(MS) that is customary in CHPT. The precise procedure has been discussed in 



great detail in Ref. [29 



The procedure used in Ref. [|] corresponds to subtracting only the Aq terms 
present in all the integrals, including those in Ai and A2. 

As mentioned in |^ a Taylor expansion of the coefficients introduces in 



principle new parameters 5i via the Laurent-expansion of the Li = a^/ e + Lf'^^ + 
e6i + ■ ■ ■. We have checked that the terms involving 6i take the form of a local 
action for the quantities considered in this manuscript, thus they can be absorbed 
in the LECs as proven in general in . 



We have defined 

-1 



L^ ^ (/ic)-^^ ( :^T, + L[(/.) ) = (/i)-2^ ( T^T^r.Ao + L[(/.) + 0{e) ) . (126) 



In the main text we have suppressed the explicit /i-dependence of the L[. The 
coefficients are given in Ref. Ill] and Inc = -l/2(ln(47r) - 7 + 1). The order 



e term in the last part of Eq. ( p.26| ) has been used as well to check the explicit 
cancellations of ^(4-7?) and 7 in all expressions. 

Similarly the coefficients in the Lagrangian are used to absorb the remaining 
infinities via 



Q = (/ic)-^^ + ^ + CKf^)) = /^"'^ il2iX2 + luXi + C[(/i) + 0{e)) . 

(127) 

Dropping the terms with Aq, Ai, Aq, replacing the Cj by C[ in the main text 
and subtracting the terms proportional to C, Aq, Ai and A2 in the expressions 
for the integrals given in the preceding appendices, gives the results in the MS 
scheme. 
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